Immunization of Susceptible-Infected Model on Scale-Free networks 



O 
O 

(N 

■4— » 

o 
O 

oo 



^ ■ 

Oh' 
i ' 

o 



CO 

o 

CO 

Oh: 



> 

OO 



: 

o . 

CO ' 

o : 
• i— i 

CO ' 

^ : 

Ph. 



X 



Wen-Jie Bai, Tao ZhouE and Bing-Hong Wang 
Department of Modern Physics, University of Science and Technology of China, Hefei 230026, PR China 

(Dated: February 2, 2008) 

In this paper, we investigate two major immunization strategies, random immunization and tar- 
geted immunization, of the susceptible-infected (SI) Model on the Barabasi- Albert (BA) networks. 
For the heterogenous structure, the random strategy is quite ineffective if the vaccinated proportion 
is quite small, while the targeted one which prefers to vaccinate the individuals with the largest 
degree can sharply depress the epidemic spreading even only a small amount of population are 
vaccinated. The analytical solution is also obtained,which can capture the trend of velocity change 
versus the amount of vaccinated population. 
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I. INTRODUCTION 

Epidemic dynamics, one of the attracting problems in 
both biological and physical communities, aims at ex- 
plaining the dynamical processes of disease spreading, 
computer virus prevailing, and so on. The previous in- 
vestigations based on the differential equations are al- 
ways under the assumption of both homogeneous infec- 
tivity and homogeneous connectivity of each individual 
Therefore, if denote the possible contacts, along 
which the infection spreads, by edges, then the previ- 
ous studies mainly concentrated on the epidemic dynam- 
ics on complete or random networks. However, against 
the above assumption, the empirical data of real net- 
works indicate the universal existence of heterogeneous 
topologies y>yil3- O ne intriguing finding is that many 
real networks have approximately power-law degree dis- 
tributions, that is to say, the probability distribution 
function of degree, P(k), approximately obeys the form 
P(k) ~ fc" 7 , with 2 < 7 < 3. This kind of distribu- 
tion implies an unexpected abundance of vertices with 
very large degrees, i.e., the so-called "hubs" or "super- 
spreaders". This series of networks, named scale-free 
(SF) networks, attract many researchers to investigate 
the corresponding epidemic behaviors (see the review pa- 
per and the references therein). 

The most exemplary models in this field, having 
been extensively lucubrated, are the susceptible-infected- 
susceptible (SIS) and susceptible-infected-removed (SIR) 
models. The recent works about SIS 0,0 and SIR [alTo| 
models on SF networks present us with completely new 
epidemic propagation scenarios that a highly heteroge- 
neous structure will lead to the absence of any epidemic 
threshold. However, many real epidemic processes can- 
not be properly described by the above two models. In 
this paper, we especially focus on the onset dynamics of 
epidemic outbreaks, whose behavior is governed by the 
pure prevalence without the natural recovery or remov- 
ing. That is to say, a speeding time-scale is much smaller 
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than the recovery time-scale. In the real world, when 
the speed of the disease is so drastic that the effect of 
recovery and death can be ignored, it is the identical pro- 
cess with above one. In addition, in some broadcasting 
processes, each node is in the possession of two discrete 
state: received or unreceived. Different from the SIS or 
SIR model, the ones having received signals will not alter 
back to the state unreceived. Hence, it is more proper 
to utilize the so-called susceptible-infected (SI) model to 
describe those dynamic processes, in which the infected 
nodes stay infected and spread the infection to the sus- 
ceptible neighbors with rate A. 

Very recently, Barthelemy et al. [Til IT^ studied the 
SI model in Barabasi- Albert (BA) networks and 
found that this epidemic process has an infinite spread- 
ing velocity in the limit of infinite population. Following 
a similar process on random Apollonian networks and 
weighted scale-free networks, Zhou et al. investigated 
the effects of clustering ^ij and weight distribution [l^ 
on SI epidemics. By using the theory of branching pro- 
cesses, Vazquez obtained a more accurate solution about 
the time behavior of SI model [lfj- The SI model with 
identical infectivity, which leads to a slower spreading in 
SF networks than the standard model, has recently been 
investigated by Zhou et al. 01- And the geographi- 
cal effect on SI model is studied by Xu et al. [T^, [l!j . 
Although these previous works are very helpful to the 
deeply understanding of SI epidemic, compared with the 
extensively studied SIR and SIS models, the SI model has 
not been carefully investigated thus far. Especially, the 
immunization effect on SI dynamics, which is very im- 
portant for controlling the prevalence, has not yet been 
investigated. In this paper, we focus on the immunization 
effect of SI model on SF networks, which can be consid- 
ered as a complementary work of the previous studies on 
the immunization of SIR and SIS models. 



II. THE MODEL 

In the standard network SI model, each individual is 
represented by a node of the network and the edges are 
the connections between individuals along which the in- 
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FIG. 1: (color online) The infected density i(t) vs time. 
The four plots exhibit the time evolution of i(t) for different 
spreading rates A in normal (a) and single-log (b) coordinates, 
and for different minimal degree m in normal (c) and single- 
log (d) coordinates, respectively. The numerical simulations 
are implemented based on the BA network of size N = 5000. 
In plots (a) and (b), the average degree of BA networks is 
fixed as (k) = 6, and in the plots (c) and (d), the spreading 
rate is fixed as A = 0.001. The legends in panels (a) and (b) 
denote the different spreading rates, and the legends in panels 
(c) and (d) denote the different minimal degrees. All the data 
are averaged over 1000 independent runs. 



fcction may spread. Each individual can be in two dis- 
crete states, either susceptible or infected. The infection 
transmission is defined by the spreading rate A at which 
each susceptible individual acquires the infection from an 
infected neighbor during one time step. 

Using the mean-field the ory , the reaction rate equa- 
tions can be written as [lU 03 : 



dikjt) 
dt 



X(k)(l-i k (t))Q k (t) 



(1) 



where ik(t) denotes the density of infected individuals 
with degree k, (k) the average degree, and 9^ the den- 
sity of the infected neighbors of a fc-degree node. Ne- 
glecting terms of order 0(i 2 ), the evolution behavior, 



m = Ekik(t)p(k), 
miii : 

t(t) ~ e ct 



can be approximately solved as 



with c ~ (k 2 )/(k). 



(2) 



In an SF networks with a degree distribution exponent 
2 < 7 < 3, the second-order moment (k 2 ) will approach 
to infinite as the increase of network size N, indicating 
an infinite velocity in the onset of epidemic spreading. 

In Fig. 1, we report the simulation results about the 
time evolution of infected density with initially one ran- 
domly selected node to be infected. All the simulations 
are implemented on the BA networks [l3| . which can 
be constructed by continuously adding one node with m 



edges connected to the existing nodes relying on the prob- 
ability proportional to their degrees. The advantage with 
the BA model is that it is the mostly studied and lacks 
structural-biases such as none-zero degree-degree corre- 
lations. Clearly, the epidemic spreading if very fast, and 
in the early stage, i(t) follows an exponential form. 



III. IMMUNIZATION EFFECT 

Immunity is a practical controlling strategy to the 
prevalence of the disease. The most extensively investi- 
gated approaches is the so-called mass vaccination [Hl20| 
(or called random immunization) . In random immuniza- 
tion, a fraction / of the whole population is randomly 
selected to be vaccinated in advance. The most signifi- 
cant problem is that whether it is effective for the highly 
heterogeneous networks? In the previous works, by using 
the mean-field theory and branch process theory, Callway 
et al. 0| and Cohen et al. [H , separately but almost at 
the same time, both proved that the random immuniza- 
tion is of less effectivity for SIR model on SF networks. 

In Fig. 2, we plot the time evolution of infected density 
i{t) for different immunization range /, which is defined 
as the fraction of population being selected to be vac- 
cinated. From Fig. 2, one can find that the spreading 
velocity has almost no change if only a very few individ- 
uals are selected to be vaccinated. Therefore, similar to 
the situations for SIR model, the random immunization 
is of less effectivity for SI model on SF networks. 

Other than the random immunization, if the degree of 
each node is known, one recently proposed efficient im- 
munization strategy is the so-called targeted immuniza- 
tion which means to vaccinate the nodes with the 
largest degrees first. Fig. 3 shows the effect of targeted 
immunization for different /. The spreading velocity re- 
markably decreases even only a small fraction, / = 0.001, 
of population get vaccinated, which strongly indicates the 
efficiency of the targeted immunization. From Fig. 3(b), 
it is observed that the time scale governing the epidemic 
behavior in the early stage sharply changes even only 
10~ 3 fraction of population (i.e. five nodes) are vacci- 
nated. 

Consider a scale-free network of size N, the degree dis- 
tribution, P(k) = Ak~ 7 , obeys the following normalized 
condition 

Ak^'dk = 1, (3) 



where A is a normalized constant, M the maximal de- 
gree and m the minimal degree. We assume after the 
fraction / of population with largest degrees having been 
vaccinated, the maximal degree decreases to k c (f), and 
if / is sufficiently small so that the degree distribution 
still obeys a power-law form with exponent 7 almost un- 
changed, then 

ffecCf) 

Ak^dk = 1 - /. (4) 
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FIG. 2: (color online) The infected density i(t) vs time under 
random immunization for the SI model based on BA networks 
in normal (a) and single-log (b) plots. The network size N = 
5000, the minimal degree m — 3, and the spreading rate A = 
0.001 are fixed. The black, red, green, blue and sky-blue 
curves, from top to bottom, represent the cases of / = 0, 
0.001, 0.002, 0.004 and 0.008, respectively. All the data are 
averaged over 1000 independent runs. 



FIG. 3: (color online) The infected density i(t) vs time under 
targeted immunization for the SI model based on BA networks 
in normal (a) and single-log (b) plots. The network size N = 
5000, the minimal degree m — 3, and the spreading rate A = 
0.001 are fixed. The black, red, green, blue and sky-blue 
curves, from top to bottom, represent the cases of / = 0, 
0.001, 0.002, 0.004 and 0.008, respectively. All the data are 
averaged over 1000 independent runs. 



Following the mean- field theory [lj 

Il2| . the time evo- 
lution of i(t) in the early stage approximately obeys an 
exponential form i(t) — i(0)e T ^' t , where i(0) = 1/N is 
the initial infected density. The time scale r(/) can be 
obtained as 



r(f) = A 



- 1 



(5) 



where P'(k) is the degree distribution of the network af- 
ter vaccination (i.e. after the removal of N f nodes of 
largest degrees), which reads 



P'(k) 



1-f 



P(k), fc = m,m +!,-•• ,k c (f). (6) 



In the large- limit of N, the maximal degree in the original 
network, M ~ \A/V, approaches to infinite. Combine Eqs. 
(2), (3), (4) and (5), the time-scale r(/) after targeted 
immunization in the N — > oo limit can be analytically 
obtained, for any 7 € (2, 3), as 



r(/) = A m x 



7 /(3-7)/(l-7) 



1 



3-7 /(2-7)/(l-7) - 1 



1 



(7) 



Especially for the BA networks with 7 = 3, the analytical 
result is 



r(f) = A 



mlnf 

2(V7 : i) 



- 1 



(8) 



In Fig. 4, we report the numerical and analytical re- 
sults for BA networks with N = 10000 and m = 3. Al- 
though the analytical result for BA networks (e.g. Eq. 
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FIG. 4: (color online) The rescaled time scale t/A vs. im- 
munization fraction /. The black squares represent the sim- 
ulation result based on BA networks. The network size 
N — 10000, the minimal degree m = 3, and the spreading 
rate A = 0.001 are fixed. The blue and red curves denote the 
analytical results for 7 = 3 and 7 = 2.85, respectively. All 
the data are averaged over 1000 independent runs. 

(8), shown as the blue curve) can capture the trend of 
t(/), the quantitive departure is very obvious. Note that, 
the fitting value of 7 in finite size BA networks is smaller 
than 3.0, which will lead to a even broader distribution 
than P(k) ~ fc~ 3 thus a faster spreading than the theo- 
retical prediction. We have obtained the average fitting 
value of 7 as 7 ~ 2.85, over 100 independent configu- 
rations of BA networks with N = 10000 and m = 3. 
The red curve in Fig. 4 represent the analytical result 
for the modified exponent 2.85 following Eq. (7), one 
can see clearly that it agrees well with the simulation for 
small / and can capture the trend of r(/). For larger /, 



the assumption that the degree distribution still obeys a 
power-law form with same exponent after the removal of 
Nf hub nodes will not be valid, resulting in the observed 
departure. 



IV. CONCLUSION 

As an important branch of the studies on epidemic 
spreading, immunity never loses its attraction. Some 
striking conclusion somewhat changes our opinions about 
epidemic. However, despite of the well-studied SIS and 
SIR model, the immunization effect on the outbreaks of 
epidemic spreading, of significantly practical value, has 
not been carefully investigated thus far. The purpose of 
this paper is to provide a complementary work of the 
previous studies on the immunization of SIR and SIS 
models. 

Two major immunity strategies are investigated based 
on the BA networks. The random immunization is of less 
cffectivity while the targeted immunization can sharply 
depress the spreading velocity even only a very few hun 
nodes are vaccinated. Furthermore, the analytical results 
is obtained which agree with the simulation well for suf- 
ficiently small immunization fraction. 
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